EN FR
EN FR


Section: New Results

White Noise-based Stochastic Calculus with respect to Multifractional Brownian Motion

Participants : Joachim Lebovits, Jacques Lévy Véhel.

The purpose of this work is to build a stochastic calculus with respect to (mBm) with a view to applications in finance and particularly to stochastic volatility models. We use an approach based on white noise theory.

White Noise-based Stochastic Calculus with respect to multifractional Brownian motion

The following results may be found in [28] . Integration with respect to mBm requires stochastic spaces in which we can differentiate or integrate stochastic processes. Considering the probability space (𝕊 ' (),𝔹(𝕊 ' ()),μ) where μ is probability measure given by Böchner Minlos theorem, one can build to spaces, noted (𝒮) and (𝒮 * ) which will play an analogous role to the spaces 𝕊() and 𝕊 ' () for tempered distributions. We recall that 𝕊() is the Schwartz space of rapidly decreasing functions which are infinitely differentiable and 𝕊 ' () is the space of tempered distributions. Let us moreover note (L 2 ) the space of random variables defined on the probability space (𝕊 ' (),𝔹(𝕊 ' ()),μ) which admit a second order moment. The mBm B (h) has the following Wiener-Itô chaos decomposition in (L 2 ):

B (h) (t)= k=0 + <1 [0;t] ,M h(t) (e k )> L 2 () <.,M h(t) (d k (t) ) =e k >= k=0 + 0 t M h(t) (e k )(s)ds<.,e k >(6)

where (e k ) k denotes the family of Hermite functions, defined for every integer k in , by e k (x):=π -1/4 (2 k k!) -1/2 e -x 2 /2 h k (x) and where (h k ) k is the family of Hermite polynomial, defined for every integer k in , by h k (x):=(-1) k e x 2 d k dx k (e -x 2 ). Note moreover that M H is an operator from 𝕊() to L 2 () for every real H in (0,1) and <.,e k > is a centered random Gaussian variable with variance equal to 1 for all k in . We can now define a process, noted W (h) , from to (𝒮 * ), which is the derivative of B (h) in sense of (𝒮 * ) by

W (h) (t)= k=0 + [d dt 0 t M h(t) (e k ) (s) d s]<.,e k >.(7)

Hence we define integral with respect to mBm of any process Φ:(𝒮 * ) as being the element of (𝒮 * ) given by:

Φ(s,ω)dB (h) (s)= Φ(s)W (h) (s)ds(ω),(8)

where denotes the Wick product on (𝒮 * ). It is then possible to get Itô formulas and Tanaka formula such as

0 T f x(t,B (h) (t))dB (h) (t)=f(T,B (h) (T))-f(0,0)- 0 T f t(t,B (h) (t))dt-1 2 0 T d dt[R h (t,t)] 2 f x 2 (t,B (h) (t))dt.(9)

for functions with sub exponential growth and where the last equality holds in L 2 .

Once this stochastic calculus with respect to mBm is defined, we can solve differential equations arising in mathematical finance.

Multifractional stochastic volatility

Multifractional stochastic volatility

The results of this part may be found in [6] . We assume that, under the risk-neutral measure, the forward price of a risky asset is the solution of the S.D.E.

dF t =F t σ t dW t ,dln(σ t )=θμ-ln(σ t )dt+γ h d B t h +γ σ dW t σ ,σ 0 >0,θ>0,(10)

where W and W σ are two standard Brownian motions and B h is a multifractional Brownian motion independent of W and W σ with functional parameter h, which is assumed to be continuously differentiable. We assume that W is decomposed into ρdW t σ +1-ρ 2 dW t F , where W F is a Brownian motion independent of W σ . Note that d B t h denotes differentiation in the sense of white Noise theory. The solution of the volatility process (σ t ) t[0,T] is

σ t = a.s.expln (σ 0 ) e -θt + μ 1-e -θt + γ σ 0 t e θ(s-t) d W s σ + γ h e -θt I t B h ,(11)

where I t B h := 𝑎.𝑠e θt B t h -θ 0 t e θs B s h ds.

Since the solution the previous S.D.E. is not explicit for (F t ) t[0,T] we use preconditioning and then cubature methods in order to get an approximation of it. This model allows to take into account the well-known "smile" effect of volatility, as well as its evolution at various maturities.

Approximation of mBm by fBms

In [18] , we establish that a sequence of well-chosen lumped fractional Brownian motions converges in law to a multifractional Brownian motion. This allows to define stochastic integrals with respect to mBm by "transporting" corresponding stochastic integrals with respect to fBm.